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9? ' Abstract 



We use VIPT in calculating the elastic form factors and charge radii of D, Dg, B, Bg 
and Be mesons in a QCD inspired potential model. For that we use linear cum Coulombic 
potential and opt firstly the Coulombic part as parent and then linear part as parent. The 
results show that charge radii and form factors are quite small for Coulombic parent case 
than the linear parent. Also, the analysis leads to a lower as well as upper bound on the four 
momentum transfer hinting at a workable range of within this approach which may be 
, useful in future for experimental analysis. Comparision of both the options prefer the linear 

parent. 
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1 Introduction 

The potential model description in the nonperturbative regime of QCD is tremendously successful 
in providing both qualitative and quantitative description of the hadron spectrum and the deacy 
modes. On this basis, we have pursued a Non Relativisic Constituent Quark Model (NRQM)[l-4] 
for mesons containing a heavy (light) quark and a light (heavy) antiquark .Although non rela- 
^ ' tivistic in nature ,relativistic effect is to be introduced from outside [5,6] due to the light quarks 
involved. Basically, the model relies on the work of Rujula,Georgi and Glashow [7] who used a non- 
relativistic treatment of potential model which was quite successful in the description of different 
hadronic properties. We have solved the Schrodinger equation for the spin independent Fermi-Breit 
Hamiltonian (used in the work of Rujula et al [7]) consisting of the linear cum Coulombic potential 
for the ground state [2]. The solution i.e. wavefunction has been obtained using different approxi- 
mation methods like Dalgarno method [8] and Variationally Improved Perturbation Theory(VIPT) 
[9-11] which is then used in predicting the Isgur-Wise (I-W)function [2, 3,4, 12, 13], mass, decay con- 
stant, charge radii[l] etc. We note that with the linear cum Coulomb potential of QCD we have two 
options in choosing parent or child(i.e.perturbation)-(i)first, we can consider the linear one as the 
perturbation (i.e. the Coulombic one as the parent ) and then (ii) Coulombic one as the perturba- 
tion (i.e. linear one as the parent). As we have already successfully used VIPT for both the options 
in the calculation of slope and curvature of I-W function for D, D^, B. Bg and B^ mesons [12, 13], so 
extending it for the prediction of elastic form factors and charge radii (which provide important 
insight on the distribution of different charge constituents of a hadron) definitely makes a sense . 
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It is well known that the form factor and charge radius are dependent on the momentum 
transform of the wavefunction.So,getting an appropriate wavefunction is very essential for a fruit- 
ful analysis. With the success of VIPT in the calculation of Isgur-Wise function as pointed above 
[12,13], one can expect a similar success here also .It is worthwhile to note that while investigat- 
ing the form factor ,one must take into account of a proper range of four momentum transfer 
Q^.The range usually determines the applicability of perturbative QCD(pQCD) or nonper- 
turbative(npQCD).So,an accurate selection of range within the nonperturbative approach is 
necessary which will also fall within the experimental regime. This fecilitates a direct comparision 
between theory and experiment. This has been done both theoretically and experimentally since 
long [14-19] for the light tt, K etc mesons. However,for the mesons which contain atleast one heavy 
quark,very little have been investigated theoretically [20-22]. In the absence of any experimental 
data for them,our results may be helpful in future in the experimental set up regarding the 
range . 

As far as our model is concerned, the perturbative or nonperturbative regime of QCD can be 
interpreted [1] through the relativistic factor 'e = 1 — ^1 — ^|^' . The reality constraint on the 
form factor F (Q^) leads to the condition < e < 1, where the case e — t- ( e — t- 1) corresponds 
to the perturbative (nonperturbative) limit of QCD .The e — ?■ 1 limit demands large or low 
Q^.So, discussing the nonperturbative effects of QCD with large confinement parameter b ,we 
must consider the low limit of in this model .However, we have observed in Ref[l] that 
large value of b{= 0.183GeV^) prohibits the use of low compelling one to involve with small as 
which corresponds to the perturbative regime and thus can't be accounted in this nonperturbative 
approach. 

We reanalyze all these observations in this approach of VIPT for both the cases -linear or con- 
finement part as perturbation and Coulombic part as perturbation. We will explore the possibility 
of incorporating significant value of as even with large confinement (6 = 0.183GeV^). This work 
will also check the status of both confinement and Coulombic part as perurbation and observed the 
consequences regarding the usable range of to work with in the absence of experimental data 
for the said mesons .The calculations are done with a fixed value of as from V-scheme [23-25] with 
large confinement effect b = O.lSSGeV"^ instead of variation in both. Even with this single value of 
as and b one can draw similar conclusion regarding the effective range of Q^. The calculated form 
factors are plotted graphically to show their variation with for both the cases. 

Basically, this work explores the possibility of improving the results for form factors and charge 
radii over those of Ref[l,22] with the help of VIPT. In the process, we also try to find a useful 
range of which may be workable for the experimental investigations in the later course of 
time. Comparision of both the options is being made to arrive at a conclusion in using VIPT. 

The rest of the paper is organized as follows : section 2 contains the formalism, section 3 the 
result and calculation while section 4 includes the discussion and conclusion. 



2 Formalism 

2.1 VIPT with Coulombic potential as parent 
2.1.1 Wavefunction 

We breifiy reformulate the VIPT with the expression for the wavefunction ^lUl [12] corrected upto 
the first order of j*'* state given by : 



^3 ^ 2-^ ^,(0) i:,{0) 
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where P' is the variational parameter (which is later optimized w.r.t. energy) considered instead 
of physical parameter P. For the Coulombic part as parent,the physical parameter is a = ^ and 
the optimized variational parameter is a' [12] .As done in Ref[12],we consider the wave function for 
triple term consideration of the summation given by equation (1) above and rewrite the equation 
(viz. eq(45) of Ref[12]) with the subscript 10 (n = 1, / = 0) being replaced by T : 



^T = ^i?-A(l-^) e-*+5(l-^ + ^^)e-* + 



U 16 32 8 X 96 ' ^ ' 



where the different parameters are given by: 



c\ = ^ (3) 



A= r ^A^Q^io (q^ - «io) 326 

30F(a;o)^ 27 81Mo^ 



^_ r S/^g^o (a - g^o) 27b 

64 256^0^ 

and ^ 

^/ _ (Mo)^ r 36 (a - a'lQ) 384b 

v/? ^ 15625a;o 78125;u2a;V 
The wavefunctions for single and double term consideration can be obtained by putting B = D' = 
and D' = respectively in the same equation (2). 

However we will consider the relativistic version (e 7^ 0) of the above wave function viz. ^12j 

iJT,Rei (r) = ipT {rfia'^Qy' (7) 
The relativistic factor e is given by [1]: 



1-, i-i^ 



2.1.2 The elastic charge form factor and charge radii 

The form factor can be expressed as |26j : 

, 1 g. r+00 

=E7f / r\ijT,Rei{r)\'smQ,rdr (9) 

^ 

where 

Q. = (10) 

Putting (2) and (7) in (9) we get the form factor as : 

eF (q2) = (^iN^TiS - 2e) [q^ + q2 + + + q5 + + Qr + Qs + q9 + Qio] (H) 

where A''3^ is the same normalization constant as appeared in equation(54)of Ref[12] and the 
different qi{Qi) s(i = 1, 2, 10) are defined in the Appendix. 
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The charge radius is derived as [T] 



< > 



d{eFm) I 

= A^3^r(3 - 2e) [ri + r2 + rs + r4 + rs + re + + 
+ + rio] 



(12) 



(13) 



where the different rj s {i 



1, 2, 10) are defined in the Appendix. 



2.1.3 Status of linear potential as perturbation 



The momentum transform of equation (7) is [271 EH] 



^ eM -r(3 - 2e) [pi - P2 + + P4] 



^ TTV / '^^T,Rei {r) sinQivdr 




(14) 



(15) 



The Pi s which depend on Ql,e etc are given in the Appendix. 

If linear potential is treated as perturbation then from equation(15) above the following in- 
equality must be preserved: 



This inequality leads to a lower limit of namely Qq [Ij above which one have to use the values 
of Q^.The Qq is determined from the condition: 



Due to the quark mass dependence ,(5o s have different values and they are shown in table 3. In the 
Dalgarno method approach [1], the lower limits Qg were large and the formalism failed to account 
for large confinement effect {h = 0.183(jey^)in the nonperturbative QCD regime where values 
were taken to be large. Only in the limit 6 — )■ ,the Qq values were lowered and the formalism 
worked for low range [I]. In this method of VIPT,the values of Qq are shown to be quite small 
even with large confinement effect (6 = O.lSSGeV'^) enabling us to work in the nonperturbative 
QCD regime with large a^- 

We also note that for single term consideration only p2 exists on the RHS of inequality and 
for double term both p2 and p^ exist. We have also recorded the values of Ql for single and double 
term consideration in table 3. 

2.2 VIPT with linear potential as parent 
2.2.1 Wavefunction 

As pointed in Ref[10,13] linear parent gives rise to Airy functions. The physical parameter is b and 
the optimized variational parameter is b' .We reproduce the analogous wavefunction in this case 
also for three term consideration of eq.(l) as was for the Coulombic parent: 



Pi> P2-P3- Pi 



(16) 



Pi=P2-P3- Pa 



(17) 



4 



(P02 - Poi) ^ 

(2^1) 
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(^(b -b') <r >2,i -a < ^ >2,i) '020 (r) + 



(P03 - Poi) 

(2/.) 



-|- (^(6 - 6') < r >3,i -a<- >3,ij Vso (r) 



-|- (^(6 - 6') < r >4,i -« < - >4,i ) i>4o (r)] 



(18) 



(P04 - Poi) & 

where N" is the normahzation constant as appeared in equation(25) of Ref[13].We note that for 
single (double) term consideration of equation(l) the third and fourth term ( fourth term) is 
dropped from the equation (18) and normalization constants also changes to different one (Eq.l7 

and 21 of Rcf [13]). 

For this case also we take the relativistic version of the above wavefunction : 



The zeros of the Airy function pon is given by eq.(ll) of Ref[13] as: 

"37r(4n- 1)" ' 



Pon 



and 



r+oo 

< >n,n'^ NnNn' / 

Jo 



Ai ((2/x6')3r - pon) Ai ((2//6')3r - pon') dr 



(19) 



(20) 



(21) 



where Nn, Nn' are the normalization constants for n and n' states respectively etc. 

Like the expressions we have adopted the same values of b,b ,pon from Ref[13](see tables 1,2 
there) . 

2.2.2 The elastic charge form factor and charge radii 

Putting (18) and (19) in the definition of form factor (9) given above, the form factor is found to 
be : 



(22) 



The coefficients C, C's are given in table 4.They are of course different for single , double or more 
than two term consideration. Numerical integrations are done in getting the above result. 

The corresponding charge radius is obtained by using equation (18) and (19) in(12) which are 
recorded in table 4. 

2.2.3 Status of Coulombic potential as perturbation 

The momentum transform of (19) is: 



, . g . / 2 r+oo 

'^T,Rei\Q) = Zl^y-y^ riljT,Rei (r) sinQirdr 



(23) 



Q 



' (K+K + Pt + Ps)] 



(24) 



The p[ s {i — 1,2, ..,8) are given in the Appendix. 

If Coulombic potential is treated as perturbation then from equation(24) above the following 
inequality must be preserved: 



02 

P'l+P2+P3+P4 > Y (P'5+Pe+P7+P's) 



(25) 



This inequality leads to a upper limit of namely Ql below which one have to use the values 
of Q^.The Qq is determined from the condition: 



Q2 

P'l +P'2+P3+P'a=Y + ^6 + P7 + P's) 



(26) 



The different values of upper limit Qq ^ shown in table 5. The corresponding values for single 
and double term consideration are also shown. 




Figure 1: Variation of eF{Q ) vs Q for D ,Dg and Bc-meson with Coulombic parent 



3 Calculation and Results 

We have listed b' in table 1 while the lower and upper limit of Ql for single,double and triple term 
consideration are given in tables 3 and 5. In table 2 , we record the charge radii for single,doublc 
and triple term consideration for Coulombic potential as parent; whereas the same is recorded 
for linear potential as parent in table 4. The infinite mass consideration shown by the subscript 
00 is also included for triple (single) term consideration for Coulombic (linear) parent. The table 
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I ' ' ' ' ' 

0.6 1 1.5 2 2.5 3 

Figure 2: Variation o^eF{Q^) vs for D ,Ds and i?c-nieson with linear parent . 

6 shows charge radii of different mesons obtained from other models and data. The values 
are taken from the y-scheme [221 [251 125] and the integrations are done numerically for all these 
calculations. The graphs show the variation eF{Q'^) vs for D, Dg and mesons for both the 
options. 
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Table 1: Values of 6' with b = 0.183GeV^ for relativistic case only.The values are the same as 
recorded in table 2 of Ref[13] 



Mesons 


Reduced mass ^ 


a — ^ 6' with relativistic effect 


D 


0.2761 


0.924 


16.24 


Ds 


0.368248 


0.924 


19.8 


B 


0.31464 


0.348 


5.587 


Bs 


0.4401 


0.348 


5.954 


Be 


1.1803 


0.348 


8.103 



Tabic 2: Values of charge radii for different mesons with Coulombic parent for single, double 
and triple terms in eq.(l).The subscripts 'S,D,T^ correspond to single, double and triple terms 
respectively whereas 'F' means finite mass consideration.The infinite mass limit ( subscript oo is 
used) is shown for the triple term alone . 



Meson 




D+ 


Dt 


B+ 


B° 


B'. 


Bt 


< r| p >^ infm 


-0.121 


0.115 


0.11 


0.2545 


-0.1822 


-0.168 


0.108 


< rjjp >2 infm 


-0.119 


0.112 


0.106 


0.2512 


-0.1788 


-0.164 


0.105 


< r^p >2 infm 


-0.118 


0.109 


0.101 


0.2464 


-0.1736 


-0.158 


0.1034 


< 4,oo >^ i^f^ 


-0.131 


0.12 


0.113 


0.263 


-0.186 


-0.1742 


0.1325 



Table 3: Values of lower limit of four momentum transfer Qq with Coulombic parent taking 
single, double and triple terms in eq.(l).We have to use values above these. 



Meson 


D+ 


D- 


Dt 


B+ 


5° 


B's 


Bt 


Qo,s 
Qo,D 
Qo,T 


0.0004 
0.00036 
0.0003 


0.0004 
0.00036 
0.0003 


0.001 
0.0009 
0.0007 


0.053 
0.052 
0.05 


0.053 
0.052 
0.05 


0.075 
0.072 
0.07 


0.211 
0.209 
0.205 
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Table 4: Values of charge radii for different mesons with linear parent for single, double and triple 
terms in eq.(l) . The subscripts '5", D, T" correspond to single,double and triple terms respectively 
whereas 'F' means finite mass consideration .The infinite mass limit ( subscript oo is used) is 
shown for the single term alone. 



Meson 




D+ 


Dt 


5+ 


5° 






Cs 


8.24 


8.24 


5.916 


14.22 


14.22 


10.91 


4.50 


C's^ 


2.266 


2.266 


1.173 


7.71 


7.71 


4.52 


0.78 


< rgp >2 infm 


-0.197 


0.1494 


0.104 


0.425 


-0.2996 


-0.2227 


0.1125 


' Cd 


13.1 


13.1 


8.1 


26.7 


26.7 


16.5 


5.2 




2.69 


2.69 


1.67 


9.89 


9.89 


5.7 


1.1 


< r^ p >2 injm 


-0.21 


0.161 


0.121 


0.473 


-0.336 


-0.2489 


0.127 




18.14 


18.14 


11.18 


43.09 


43.09 


25.56 


6.636 


Crp 


3.365 


3.365 


2.199 


13.14 


13.14 


7.655 


1.363 


9 1 • C 

< r^ p >2 injm 


-0.24 


0.182 


0.143 


0.555 


-0.391 


-0.289 


0.148 




-0.246 


0.174 


0.125 


0.453 


-0.32 


-0.246 


0.144 



Table 5: Values of upper limit of four momentum transfer Qq with linear parent taking sin- 
gle, double and triple terms in eq.(l).We have to use values lower than these. 



Meson 


D+ 


D- 


Dt 


B+ 




B's 


Bt 


^0,5 

Qo,D 
Qo,T 


2.297 
2.1 
1.88 


2.297 
2.1 
1.88 


2.92 
2.67 
2.387 


1.43 
1.31 
1.177 


1.43 
1.31 
1.177 


1.676 
1.56 
1.386 


3.11 
2.89 
2.55 



9 



Table 6: Prediction of < >2 infm for finite and infinite mass consideration in other models. Tlie 

subscript 'F' ('oo' ) means finite (infinite) mass limit. 

Meson D° D+ D+ B+/B- B^JB^^ B+/B~ 

< rl >5[9] I 0.506 0.491 0.258(5-) I 0.256(5°) 0.236(5-) 

<r|.>5[10] -0.551 0.43 0.352 0.612(5+) -0.432 -0.345(5°) 0.207(5+) 

<r^>5[10] -0.704 0.498 0.425 0.704(5+) -0.498 -0.425(5°) ...(5+) 

<4>5[11] -0.484 0.366 0.355 1.72(5+) -1.21 -1.17(5°) 1.43(5+) 

<r^>5[ll] -0.6025 0.427 0.427 1.836(5+) -1.29 -1.29(5°) 1.84(5+) 



4 Discussion and Conclusion 

Wc have analyzed elastic form factors and charge radii in a QCD inspired potential model using 
VIPT under two scenarios-Coulombic potential and linear potential as parent . 
We summarize our achievements below: 

I. The form factor eF(Q'^) decreases with the increase of (as it should) for both the scenarios. 

II. The form factor is either very small (for 5-sector mesons)or small ( for 5-sector mesons)with 
Coulombic parent as compared to those with linear parent. 

The charge radii is also observed to be smaller with Coulombic parent as compared to linear 
parent. 

III. We use a fixed set of values for ag under y-scheme[23-25] in the calculation , for example 
it is 0.693 for the D, Dg mesons which is larger than the value 0.261 for the 5, Bg, Be mesons. This 
consideration directly results in the unexpected smaller values of charge radii for D, Dg mesons as 
compared to the B, Bg, Be mesons. Larger ag values are responsible for smaller charge radius. 

IV. While checking the status of confinement as perturbation i.e. Coulombic parent (or 
Coulombic part as perturbation i.e. linear parent) we end up with a lower (or upper) limit 
on Q"^ .This allows us a useful range of Q"^ to show the variation of form factor which is shown in 
fig 1 and 2. 

V. In the present analysis, even with large b,the lower limit of Q^{foT linear perturbation) are 
really small as shown in tabic 3 for fixed a^-We have seen that for ag = 0.693 , the lower limit of 

for D,Dg mesons are respectively 0.0003,0.0007, whereas with ag = 0.261 ,the lower limit of 
for 5, Bg, Be mesons are respectively 0.05, 0.07, 0.205 .These values for 5, Bg, B^ mesons will be 
lowered if we put ag > 0.261. This is clearly advantageous over the Dalgarno method with linear 
perturbation as done in Ref[l] where the formalism broke down for large fe.Thus, this approach 
allows a large vahic of ag{Q'^) in the limit ^ even with large confinement , an important 
feature absent in Ref[l]. 

VI. The Coulombic perturbation leads to an upper limit of Q"^ in this case (table 5) which 
allows us to use any value of as(Q^) in the limit — )■ 0. 

VII. Further, if we look at Eq.(l), consideration of different terms leads to different charge 
radii and the limiting values of for both the cases. The charge radii and the lower limit of 
decrease with more terms for the linear part as perturbation (i.e. Coulombic parent) whereas the 
charge radii increase and upper limit of decreases for the Coulombic part as perturbation (i.e. 
linear parent). 

VIII. The infinite mass consideration in this work shows that the charge radii are larger than 
those for finite mass consideration to agree well with other models (table 6). 
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The above list as a whole suggests success of VIPT over the Dalgarno method [H |22] as far 
as large confinement and limiting values of are concerned. The difference in the values of form 
factors and charge radii for both the cases may be attributed to the use of same (i.e. Q^) under 
y-scheme for both the scenarios as the Coulombic potential is dominant for large (i.e. low r)and 
the linear potential in the low (i.e. large r)regime.We must note that we have used the low 

assumption in the calculation of form factors and this clearly effects the upper limit of 
corresponding to the validity of Coulombic perturbation (i.e. linear parent). The larger value of 
for D-sector as compared to S-sector is also another point to take account of this .Although, the 
linear parent has shown more flexibility and hence is the better option than Coulombic parent in 
VIPT, but it has used terms up to a particular order in'r' in the integration involved with Airy 
function (which is an infinite series). This may lead to loss of certain information as far as physics 
is concerned. In the absence of any experimental results for these mesons , it is quite difficult 
to make a direct conclusion but there is clear indication that one must be careful in choosing 
the parameter q;s(Q^) as well as the confinement parameter in the calculation of form factor and 
charge radius within the QCD framework. 

The above discussion led to the conclusion that there is scope to use this approach in the 
study of meson decays. The lower and upper limit on (i.e. range of Q^) in this analysis may be 
useful in the experimental set up to investigate cross-section and form factor in future for these 
mesons. Further, from the model specific values of form factors and charge radii the method allows 
to investigate the behaviour of w.r.t in the nonperturbative regime of QCD. 
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Expressions for s : 



Qi 



(4/i2a'^+g2)(l-e) 



1 



(3 - 2e)fla' ^{A- 2e) (3 - 2e)/i V 



(/X2a''+Qf)(l-^) (//2a'^ + gf)(l-5-^) (/X2a'' + Q2)(2-6) 



(27) 
(28) 



4(3 - 2e)^a' 



i^ + QlY'-'^ 3(^ + g2)(i.5-.) 



16(4 - 2e)(3 - 2e)^2^' 8(5 - 2e)(4 - 2e)(3 - 2e)n^a' 



27(^ + gf)(2-^) 



81(^ + gf)(2-5-e) 



+ 



+ 



4(6 - 2e)(5 - 2e)(4 - 2e)(3 - 2 e) ix^oi 
729(^ + gf)(3-) 



(29) 



12 



54 = D'' 



16(^ + Qf)(i- 



3(3 - 2e)na' 

32(^ + gf)a-5-^) 



17(4 - 2e)(3 - 2e)^Ji'^a' 19(5 - 2e)(4 - 2e)(3 - 2e)^Ji^a' 



256(^ + Qf)M 



1536(^ + gf)(2-5-^) 



+ 



+ 



7(6 - 2e)(5 - 26) (4 - 26) (3 - 2 e) 
6144(^ + g2)(3-6) 



(7 - 26) (6 - 26) (5 - 26) (4 - 2e)(3 - 2e)ii^a' 



q5 = 2c[A 



12288(^ + g|)(3-5-^) 
1 (3 - 2e)fj,a' 



qe = 2c[B 



1 



+ gf)M 2(^ + gf)(i-5-) 



2(3-26)//a' ^ 2(4 - 26)(3 - 26)/i2a'' 



q, = 2c[D'[ 



g + g2)(i-5-e) 27(^ + g?)(2-) 

3(3 - 2e)/xa' 



4(^ + gf)a- 

,2 



16( 



I /02\(1.5-e) 
16 ' -"V 16 ^'^i^ 

(4 - 2e) (3 - 2e)// V (5 - 2e) (4 - 2e) (3 - 2e)// V , 



32 ( 



16 



768 ( 



16 



g2)(2.5-e) 



+ 



(30) 
(31) 
(32) 



(33) 



-2AB[ 



1 _ 5(3 - 2€)fM' 

+ QIY'-'^ ~ eC-^ + g|)a-^-") 

20(4-2e)(3-2e)//V _ (5 - 2e)(4 - 2e)(3 - 26)//=^^^' ^ 

(27^^ + gf)(2-) 



27(^5ig^ + g2) (2-5-6) 



+ 



(34) 



= -2AD'[ ^ 5(3-26W_ 

^C-^ + QlY'-'^ i6(5^ + gf)(i-5- 

(4 - 2r)(3 - 2r)//--'' _ 13(5 - 2f)(l - 2r)(3 - 2r)// 
8(^ + Q?)M 



'TV' 



768(5^ + gf)(2-5-^) 

(6 - 2£)(5 - 2e)(4 - 2e)(3 - 26)//%^" 

1536(9i^g:l+g2)(3-e) 



+ 



(35) 



q,o = 2BD'[ 



1 



25{3-2e)fM' 

^4(^ + gf)(^-) " 48(^5gf:^ + gf)a-5"^ 

151(4 - 2e)(3 - 2e)//V 27.66(5 - 2e)(4 - 2e)(3 - 2e)//%'' 



864(^ + g|)(2-) 



144 +Qir-'-'^ 

3.66(6 - 26) (5 - 26) (4 - 26) (3 - 26) 
1152(25^ + gf)(3-) 
(7 - 26) (6 - 26) (5 - 26) (4 - 26) (3 - 2e)//%''. 
10368(25^ + g2)(3-5-^) 



+ 



+ 



(36) 
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Expressions for s: 

r, = 3c^ (1 + !^^)-2(4^V)-2 (2 - 2e) (37) 



.2/1 I "^i\-2r/o ^/l ,.2Tr^^\e-2 Q,,Tr'^Q 0,^2 / , ,2— \ e-2.5 



r2 = 3A'(1 + ^)-1(2 - 2e) (V"' )^-' - 3//cE'(3 - 2e)'(/a' )^-'-^ + 

0.75//V'(4 - 2e)2(3 - 2e)(/iV')^-3] (38) 



rs = 3B^(1 + ^)-l(^)^-^ (2 - 26) - 4/.a'(3 - 2e)^(^)-^-^ + 

-^(4- 26)2(3 -26)(-^)-^ - 
8^'"' .(5 _ 26)2(4 - 26)(3 - 26)(^^)^-3-5 + 



81 ' /V /V /V 9 
^(6 - 26)2(5 - 26)(4 - 26)(3 - 26)(^)-1 (39) 

n = SD'\1 + !?!i)-2[(/f!^)e-2 (2 _ 2e) - ^(3 - 26)2(^)^2.5 + 
TTT/j 4 32 4 

^_(4_2e)2(3-26)(^)-3 _ 



^^(5 - 26)2(4 - 26)(3 - 26)(^)-^-^ + 
21^'^" (6 - 26)2(5 - 26)(4 - 26)(3 - 26)(^)^- - 



6144 ' ' ' ''4 



12288 " ^'"^^ " ^'^^^ " ^'^^'^ " ^'^^^ " 26)(^)-^1 (40) 
r5 = 2c;A(1 + ^)-2[3(^)-2 (2 - 26) - 1.5/.cE'(3 - 26)2(^)-2.5] (41) 



rrij 



, _2r„/ Ida* CK \e- 



r6 = 2c;S(l + ^)-2[3(^^)-2(2-26) - 



9 



2^3 - 26)2(i^^)-2.5 + 

y 



g (4- 26)2(3 -26)(^^)-^] (42) 

r. = 2c;iP'(l + ^)-2[(g^)--23('-'^) - 
Tfi/ j 16 4 

9ai«'(3-26)2^25^^,_2.5 ^ 



16 ^ 16 

3u2ck' , ,25iu2a' 

^(4-2.)=(3-2e)(^) 

27T7'^ 



e-3 



^(5 - 26)2(4 - 26)(3 - 26)(^^)^-^-^] (43) 
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m. 



36 



rs = -2AB{1 + — ) 



M-2r 



\c-2 



3 (2 - 2e) 



2.5a.«'(3 - 2erC-^Y-'-' 



+ 



9 



-(4-2e)'(3-2e)( 



fx a 



(5-2e)2(4-2e)(3-2e)( 



36 
2Va^ 
36 



(44) 



m 



rg = -2AD'{1 + — )-2[0.75( 



rrii 



16 



16 



-(3-2e)^ 



ve-2.5 



16 



+ 



^ -(4-2e)^(3-2e)(- '"'^ 



8 



13^%' 
256 



-(5-2e)^(4-2e)(3-2e)( 



16 



\e-3.5 



16 



+ 



fj, a' 
512 



■(6 - 2e)2(5 - 2e)(4 - 2e)(3 - 2e)( 



16 



(45) 



no = 25D'(1 + ^r[0-7H^^r' (2 - 26) - 
1^(3-26)^(1^)— + 
^^(4-2e)^(3-2e)(4^)- - 



768 ' ' ' 144 



(6 - 2e)^(5 - 2e)(4 - 2e)(3 - 2e)(- 



\e-4 



3456 ' ' ^ 



^(7 - 26)^(6 - 2e)(5 - 26)(4 - 26)(3 - 26)(l^)-^i (46) 
Expressions for pi s: 

4 



= 4 ,.... - -S^^] (48) 



'^ + gf)^ 2(^ + gf)-? 



1 0.67(3 - 2e)/ia' 

~ Lu2s/2 (3-6) ~ , „2q'2 (4-e) 

(^ + Qf)~ (^ + Qf)~ 



2(4-2e)(3-2e)//%'' 
27(^ + gf)^ 
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(49) 



0.25 3(3 - 2e)iM' 



Expressions for p[ s : 



(4-2e)(3-2e)AtV 
32(M|1 + Qf)^ 
(5 - 26) (4 - 26) (3 - 26)// V 

768(^ + g?)-^ ^ ^''^ 



2 

\ / 


(51) 


2 

P2 = ^2 X (6^) ' 


(52) 


2 

Ps = X (feV) ' 


(53) 


2 

P4 = 77.4 X (ftV) ^ 


(54) 


2 

Ps = X (ftV) ' 


(55) 


2 

p'e ^rieX (b' fx) ' 


(56) 


2 

P7 = n7 X (^V) ' 


(57) 


2 

p's = '^s X (feV) ' 


(58) 



Each of the constants ni, n2, ng are different for different mesons and they have been obtained 
by numerical integration. 
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